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& ' Abstract 

In the first part of this paper (Sections EEJ), the main concern is 
with the boundary of the pseudospectrum of a matrix polynomial and, 
particularly, with smoothness properties of the boundary. In the sec- 
S^ . ond part (Sections loTI5|) . results are obtained concerning the number 

*-j | of connected components of pseudospectra, as well as results concern- 

ing matrix polynomials with multiple eigenvalues, or the proximity to 
such polynomials. 
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1 Introduction 

This paper falls into two parts. In the first (Sections I2HU), the main concern 
is with the boundary of the pseudospectrum of a matrix polynomial and, 
particularly, in view of its importance for boundary-tracing algorithms, with 
the smoothness properties of the boundary In the second (Sections EHHJ) , 
we further develop analysis begun by two of the present authors (see |12j ) 
on qualitative aspects of the pseudospectrum. This part is also influenced 
by earlier work on pseudospectra for standard eigenvalue problems by Alam 
and Bora in |2]. In particular, results are presented concerning the number 
of connected components of the pseudospectrum and proximity to systems 
with multiple eigenvalues. 

Let us begin with some formal definitions. First, a matrix polynomial is 
a function P : C — > C nxn (the algebra of all n x n complex matrices) of the 
form 

P(A) = P m X m + Pm-xX™- 1 + ■ ■ ■ + A A + P , (1) 

where A is a complex variable and Pq, Pi, . . . , P m G C nxn with det P m 7^ 0. 
The spectrum of such a function is cr{P) := {A G C : det P(A) = 0}. 

Since det P m 7^ 0, o~(P) consists of no more than nm distinct eigenvalues. 
A nonzero vector x G C n is known as an eigenvector of P(A) correspond- 
ing to an eigenvalue Ao G o~(P) if it satisfies P(Ao)xo = 0. The algebraic 
multiplicity of a Ao G o~(P) is the multiplicity of Ao as a zero of the scalar 
polynomial detP(A), and it is always greater than or equal to the geometric 
multiplicity of Ao, that is, the dimension of the null space of the matrix P(Ao). 
A multiple eigenvalue of P(A) is called defective if its algebraic multiplicity 
exceeds its geometric multiplicity. 

We let V m denote the linear space of nxn matrix polynomials with degree 
m or less. Using the spectral matrix norm (i.e., that norm subordinate to 
the Euclidean vector norm), we may define the max norm on V m , 

||P(A)|| := maxo<*< m ||Pi||. ( 2 ) 

Using this norm, we construct a class of matrix polynomials obtained from 
P(A) in (£Q) by perturbation. The admissible perturbations are defined in 
terms of a real polynomial w(x) = Y2T=o w j x '' with nonnegative coefficients 
and a positive constant coefficient; Wj > for each j = 1,2, ... , m, and 
Wo > 0. First consider matrix polynomials in V m of the form 

Q(X) = (P m + A m )X m + ■ • • + (P + A0A + (P + A ), (3) 



where the matrices Aj G C nxn (j = 0,1, ... ,m) are arbitrary. Then, for a 
given e > 0, the class of admissible perturbed matrix polynomials is 

B(P,s,w) := {Q{\) : IIAj-ll <swj, j = 0,l,...,m}. (4) 

This is a convex compact set in the linear space V m with the norm (J2J). 

The e-pseudospectrum ofP(X) with respect to w(x) (introduced by Tisseur 
and Higham |17|) is then 

A £ (P) := {fiEC: det Q(p) = for some Q(X) G B{P, e,w)}. (5) 

As w(x) is generally fixed throughout this paper, it will not appear explicitly 
in this notation, and we will refer to A e (P) simply as the e-pseudospectrum 
of P(X). Note that if ew m < HP^ 1 ]! -1 , then all matrix polynomials in 
B(P,e,w) have nonsingular leading coefficients, and this ensures that A £ (P) 
is bounded (Theorem 2.2 of [T2]). 

If we define the standard eigenvalue problem as that in which P(A) = 
IX — A, then it is natural to define weights W\ — (no perturbation of 
the coefficient I is admitted) and wq = 1. Thus, w(x) = 1 and, using (jSJ), 
we obtain the relatively well-understood "e-pseudospectrum of matrix A" , 
namely, 

A £ (A) = A £ (P) = {/i e C : det(//i - (A + A )) = 0, ||A || < e} . 

2 The singular value functions 

For any A G C, the singular values of a matrix polynomial P(A) are the 
nonnegative square-roots of the n eigenvalue functions of P(A)*P(A). They 
are denoted by 

si(A) > s 2 (A) > ••• > s n (X) > 0. 

The real- valued function s n : C — ► [0,oo), given by the smallest singular 
value, provides more information about the matrix polynomial P(A) than 
o~{P) alone. This will become clear in the forthcoming section when we 
discuss the pseudospectrum of P(A). Let us first describe some general prop- 
erties of s n (X). 

It is clear that an alternate definition of the spectrum of a matrix poly- 
nomial P(A) is: 

a(P) = {AGC: s n (X) = 0}. 



The connection between the zeros of s n (X) and the eigenvalues of -P(A) can 
be made more precise using the singular value decomposition. 

Proposition 1 An eigenvalue Xq G c(P) has geometric multiplicity k if 
and only if 

si(A ) > s 2 (A ) > ••• > s n _ fc (A ) > s n _ fc+ i(A ) = ••• = s n (A ) = 0. 

Our analysis depends on an important, concise characterisation of the e- 
pseudospectrum in terms of the function s n (X). This was obtained by Tisseur 
and Higham (Lemma 2.1 of [TTj). 

A £ (P) = {AGC: s n (X)<sw(\X\)}. (6) 

Clearly, a(P) = A (P) C A e (P) for any e > 0. Thus, A e (P) is nothing 
but the level set at height of the real- valued function s n (X) — ew(\X\) , 
or that at height e of the function s n (X) U7(|A|) _1 . Notice also that in the 
standard eigenvalue problem, eiy(|A|) = e in (0). More generally, eiu(|A|) 
(in equation (jUJ)) is a radially symmetric non- decreasing function of A. 
By using the Euclidean vector norm, 

• \\PW U \\ 
s n {X) = mm . . . (7) 

u^O \\u\\ 

Our first theorem has been originally established by Davies in the more gen- 
eral context of holomorphic families of bounded operators. A proof is in- 
cluded here for completeness. 

Theorem 2 Let P(A) be invertible on a domain U . Then s n (A) _1 is a sub- 
harmonic function on U . 

Proof. First we recall one of the characterisations of continuous subharmonic 
functions (see Ahlfors pQ, for example). A continuous function : U — > R 
is subharmonic if and only if, for any closed disc in U with centre Ao and 
radius r, 

r-27T 



0(A O ) < ^- [* <P{X + re ie )d6. 
2vr Jo 



A well known result from operator theory establishes that for any bounded 
linear operator T on a Hilbert space, 

Re(T0,^) 
||T|| = sup ■ 



If the Hilbert space is finite dimensional, it is easy to see that the supremum 
is attained. By virtue of (J7j), s n (X) = ||P(A) _1 || _1 . Thus, 



s n {\) 



Re(P(A)" 1 u, v) 
max Huhm 

u,v^0 \\u\\ \\v\\ 



(8) 



Now note that s n (X) l is continuous on U, and let Xq Elf and uq, Vq G C n 
be such that 



s n (Ac 



! Re(P(A ) ^o, v ) 



\\ u o\\ \\Vo\\ 
The function (P(A) _1, uo, vo) is analytic on U and so the real function 

,/n Re(P(A)- 1 Mp, Vq) 

\\ u o\\ Foil 

is harmonic on U. Furthermore, it follows from (jHJ) and © that h(X) < 
Sn(A) -1 on U. Consequently, 

Sn(Ao)- 1 = h(X ) = ^- [\(X Q + re i9 )d9 < ±- [ \ n (X + re i9 )- x d9, 

^ h 2lT Jo 

and the result follows. □ 

An important characteristic of subharmonic functions is the fact that they 
satisfy the maximum principle. Therefore, the only local minima of s n (X) are 
those A G a(P). 

The subharmonicity of s n (A) _1 has been considered recently by various 
authors. In |3j, Boyd and Desoer discuss this property in the context of 
linear control systems. Concrete applications of this theorem may be found 
in jH] for the linear case, and in [HHH] for the quadratic case. In [3], the result 
is applied in support of a certain novel procedure for finding eigenvalues of 
self-adjoint operators in infinite dimensional Hilbert spaces. 

Corollary 3 For all e > 0, every connected component of A £ (P) has non- 
empty interior. 

Proof. Suppose, on the contrary, that Q is a connected component of A £ (P) 
with empty interior. Since u>(|A|)/s n (A) is subharmonic, minA e g[s n (A)/u>(|A|)] 
should be attained at all points of Q. Thus, necessarily, Q should be a single 



point and in fact one of the eigenvalues of -P(A). The continuity of s n (X) and 
the fact that w(0) = w > ensure that this is not possible. □ 

In general, s n (X) itself is not a subharmonic function as it does not satisfy 
the maximum principle (a concrete example may be found at the end of this 
section). However, as we will see next, s n (x) is locally regular. 

First consider the nonnegative eigenvalue functions generated on C by 
the matrix function P(A)*P(A), say Si (A), 62(A), . . . , S n (X). They can be 
organised in such a way that they have a strong smoothness property. 

Lemma 4 For any given analytic curve £ : R — ► C, the eigenvalues of 
P(X)*P(X) can be arranged in such way that, for all j, Sj(((t)) are real 
analytic functions of t EM.. 

Furthermore, if s n (X) = min^S^A)) 1 / 2 is a non-zero simple singular 
value of P(X) and u\,v\ are associated left and right singular vectors, re- 
spectively, then (writing X = x + iy) s n (-) is a real analytic function in a 
neighbourhood of X and 

^ ,,., f R ( * dP{x + iy) \ ( , dP{x + iy) \\ 
Vs n (x + iy) = I Re I u x v x J , Re I u x v x I I . (10) 

The first statement follows from Theorem S6.3 of [9 j (see also Theorem 
II-6.1 of JDl)- For the second and third, see [15] . for example. 

We can interpret the first part of this lemma pictorially in the following 
manner. For £ E R, the graphs of Sj(((t)) (j = 1,2, ... ,n) are smooth and 
might cross each other. At a crossing point, the graph of the corresponding 
singular value Sk{((t)) is continuous but it changes from one smooth curve 
to another with a possible jump in the derivative (see §11-6.4 of 10J). 

We may also consider regularity properties of s n (X) as a function defined 
on the complex plane. In this case, some rudimentary ideas from algebraic 
geometry assist in discussing the n surfaces in R 3 which are (in general) 
generated by the singular values. (Where possible, the terminology of Kendig 
[TT| is followed). Write A E C in real and imaginary parts; A = x + iy, and 
define n subsets of R 3 : 

T,j := (x,y,Sj(x + iy)) ; j = 1,2, ...,n. 
Proposition 5 The union U?=i ^i ^ s a rea ^ a ^9 e ^raic variety. 



Proof. Define the function 

d{x,y,S) := det {IS- P{x + iy)*P{x + iy)) ; x,y,S eR. (11) 

Since the matrix P{x + iy)*P(x + iy) is hermitian, d(x, y, S) is a polynomial 
in x, y, S with real coefficients, and since 

n 

QSj = {(x,y,S)ER 3 : d(x,y,S) = 0}, 
i=i 

the result follows. D 

In spite of this proposition and LemmaEJ the existence of an arrangement 
of the eigenvalues of P(X)*P(X) such that the n surfaces Ej e R 3 are smooth 
everywhere is not guaranteed in general. Consider the following example. For 
the linear matrix polynomial P(X) = IX — A, where 



.4 



3/4 1 1 

5/4 1 
-3/4 



Si U E 2 U E 3 has a conic double point at (0, 0, 5/16). Therefore, no arrange- 
ment of the singular values exists ensuring Ei, E 2 and E 3 are simultaneously 
smooth at A = 0. Moreover, in this example, s 3 (0) = s 2 (0) = a/5/16, so 
note that the hypothesis of non- degeneracy of the fundamental singular value 
in the second part of Lemma H] is essential. 

For linear polynomials, the occurrence of isolated singularities in (J" =1 Ej 
is rare. In the above example the matrix A had to be carefully crafted to 
allow the conic double point around the origin. Any slight change in the 
coefficients of A would eliminate this degeneracy. 

The following useful proposition is an immediate consequence of LemmaEJ 

Proposition 6 If Sj(X) = Sk(X) for j ^ k and for all X in a non-empty 
open set O, then O = C. 

Thus, different surfaces Ej can intersect only in sets of topological dimen- 
sion at most one. 



3 The pseudospectrum and its boundary 

Now we put these ideas into the context of the study of pseudospectra. 

A fundamentally important case is that in which e > is so small that 
A e (P) consists of "small" disconnected components, each one containing a 
single (possibly multiple) eigenvalue of -P(A). As e is increased from zero, 
these components enlarge, collide and eventually intersect in various ways 
so that the boundary of A £ (P), say dA £ (P), becomes more complex. In an 
earlier paper [T2] two of the present authors studied some basic properties 
of A £ (P) and dA e (P) in support of a curve-tracing algorithm for plotting 
dA £ (P). 

Let 

F E (x,y) = F £ (x + iy) := s n (x + iy) -ew{\x + iy\) ; x,y 6 R. (12) 

Since this function is continuous in A = x + iy G C, it follows from (jUJ) that 

<9A £ (P) c {AeC: P £ (A)=0}. (13) 

Moreover, as long as s n (A) 2 is a simple non-vanishing eigenvalue of P(A)*P(A), 
differentiation in the direction of the boundary will be well-defined as a con- 
sequence of Lemma EJ However, when s n _i(A) = s n (A), this smoothness of 
the boundary may be lost. Hence our interest in the set of A G C for which 
s n (X) is multiple; curve tracing algorithms are prone to fail around these 
points, as the directional derivatives along dA e (P) may not be well-defined. 
Even though it is quite rare 1 , in general, the right side of (|13|) might 
include points in the interior of A e (P) . This can be observed as a consequence 
of either of the two unlikely events: 

(i) the surface s n (X) having a local (but not global) maximum, 
(ii) at least three multiple sheets of Uj=i ^j intersecting in a single point. 
Demmel's matrix 

A = 



-1 


-b 


-6 2 





-1 


-b 








-1 



1 This is a rather delicate point, and it seems to have been missed in the work of several 
preceding authors as in J2] and JTj. In particular, Corollary 4.3 of [2] seems to be false 
as it stands. On the other hand, this fact seems to have little, if any impact on the design 
of algorithms. 
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with b » 1 illustrates (i) for the standard eigenvalue problem with w(x) = 1. 
Indeed if b = 100 and -P(A) = (XI — A), 53(A) has a local maximum at A = 0, 

cf. m. 

Higher order examples typifying (i) can also be easily constructed. Con- 
sider, for instance, the polynomial P(A) = (A 2 — 1)(A 2 — i) in C and the 
weight function w(x) = 4x 2 + 1. The point A = is a local maximum of 

the function 

si(A) _ |A 2 -l||A 2 -z| 

w(\X\) ~ 4|A| 2 + 1 ' 

which is smooth in C\{±1, iz 1 / 2 }. This may be verified by directly com- 
puting the gradient and Hessian of this expression at A = 0. Thus, when 
e = 1 and A lies in a sufficiently small neighbourhood M of the origin we 
have si(A) < u>(|A|), so that J\f E A t (P) and <£ <9Ai(P). However, for 

6 = 1, 

P £ (0) = Fi(0) = Sl (0)-w(0) = 0, 

so in this case the inclusion of (fT5j) is proper. 

To confirm (ii), recall Example 3.5 of [2]: for w(x) = 1 and any e > 0, 
the £-pseudospectrum of P(A) = diag{A — 1, A + l, A — i, X+i} is the union of 
four closed discs with centres at the eigenvalues 1, —1, i, —i and radii equal 
to e. Thus, for e = 1, the origin lies in the set {A G C : Pi (A) = 0} but it 
is an interior point of Ai(P). 

The next result shows that dA £ (P) is made up of algebraic curves. This 
is a comforting property in the sense that the number of difficult points, such 
as cusps or self-intersections, is limited. (See Proposition 6.2.10 of |3] for an 
explicit statement of this kind.) 

Theorem 7 Let e > and assume that A £ (P) 7^ C. Then the boundary of 
A e (P) lies on an algebraic curve. In particular, dA £ (P) is a piecewise C°° 
curve, it has at most a finite number of singularities where the tangent fails 
to exist, and it intersects itself only at a finite number of points. 

Proof. We first show that <9A e (P) lies on an algebraic curve. Recall the 
function d(x,y, S) defined by (fTTj) and observe that <9A e (P) lies on the level 

set 

£1 = {x + iy : x,y G K., sw(\x + iy\) is a singular value of P(x + iy)} 

= {x + iy : x, y G M, d(x, y, e 2 w(\x + iy\) 2 ) = 0} . 



The function d(x, y, e 2 w(\x + iy\) 2 ) can be written in the form 

d(x,y,e 2 w(\x + iy\) 2 ) = a/x 2 + y 2 p(x,y) + q(x,y), 
where p(x, y) and q{x,y) are real polynomials in x,y G R. Thus, 

£1 = J a; + iy : x,|/6l, a/x 2 + y 2 p(x, y) + g(:r, y) = j . 

If u>(x) is an even function, then p(x,y) is identically zero and either £ x 
is an algebraic curve or it coincides with the complex plane. Suppose w(x) 
is not an even function. Then C\ is a subset of the level set 

£ 2 := {x + iy : x, y G K, (x 2 + y 2 )p(x, yf - q(x, yf = 0} , 

which is also an algebraic curve when it does not coincide with the complex 
plane. 

Next we show that £2 = C only if C\ = C. Thus, if £2 = C and p(x, t/), 
q(x,y) are not identically zero, then 

(x 2 + y 2 )p(x, y) 2 = q(x,y) 2 for all ij6l, 

where the order of the (irreducible) factor x 2 + y 2 in the left hand side is 
odd and the order of the same factor on the right (if any) is even. This is 
a contradiction. Hence, if £2 = C, then p(x,y) and q(x,y) are identically 
zero, and consequently, C\ = C. 

Since C\ C A e (P) and by hypothesis A e (P) 7^ C, both C±, £2 7^ C and 
so £2 is an algebraic curve. This completes the first part of the theorem. 

For the second part, note that, as s n (X) is continuous in A G C, <9A £ (P) is 
a union of continuous curves. From the above considerations it follows that 
L\ is a piecewise C°° curve and it has finitely many singularities. Then, since 

<9A e (P) c A C A e (P), 

we can actually decompose £1 = Ufc=i7fc! where 7^ (k = 1,2, ...,h) are 
suitable smooth curves with the following property: 7^ C dA £ (P) for all 
l<k<j, and 7 fc C A e (P) \ dA £ (P) for all j < k < h. Thus, <9A £ (P) = 
Ui=i 7fc as needed. □ 

Note that for the standard eigenvalue problem, w(x) = 1 is an even 
function. In this case, the above result appears in the work of Alam and 
Bora 0. 

The following technical statements will be useful subsequently. The first 
one follows immediately from (|13j) . 
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Lemma 8 If < 5 < s, then dA 5 (P) C A e (P) and A S {P) n <9A £ (P) = 0. 

In particular, note that cr(P) fi <9A e (P) = for any e > 0. 

With -P(A) as in ((TJ), consider a perturbed matrix polynomial Q(X) of the 
form (jHJ). It follows from the definition (jlj) that Q(A) G dB(P,e,w) if and 
only if ||Aj|| < ewj for each j and equality holds for at least one j. Now 
consider matrix polynomials in the interior of B(P,e,w); Int [B(P, e, to)]. It 
is easily seen that Q(A) G Int[£>(P, e, w)] if and only if 

||Aj|| < ewj whenever Wj > 0, and 

Aj = whenever Wj — 0. 

Lemma 9 // /i G dA e (P), then for any perturbation Q(X) G B(P,e,w) 
such that jjl G cr(Q), <5(A) G dB(P,e,w). 

Proof. Let /i G <9A e (P). It suffices to show that if /i G cr(Q) for a Q(X) G 
B(P, e, iy), then \\Aj\\ = ewj for some j — 0, 1, . . . ,m. Indeed, if we assume 
the converse statement, \\Aj\\ < ewj for all j, then Q(X) G B(P,e,w) for 
some e < e. But since fi G <y(Q), we have /x G Ag(P), which contradicts 
Lemma |H1 Thus, the desired assertion holds. □ 



4 The fault lines 

Differentiability along dA E (P), the boundary of the pseudospectrum, is pos- 
sible as long as the gradient of s n (X) — ew(\X\) exists and does not vanish. 
The only place where iu(|A|) might fail to have a derivative is the origin. If 
the minimal singular value, s n (X), has multiplicity one, then s n (X) is smooth 
in a neighbourhood of A. Thus, the study of those points where differentia- 
bility is lost, apart from A = 0, is confined to the region of the plane where 
the sheet of U?=i ^j corresponding to s n (X), meets the one corresponding to 
s n -i(A). This motivates the following definition. 

The rather involved indexing introduced next is required to handle mul- 
tiple eigenvalues. Below we always assume that the eigenvalues of P(A)*P(A) 
are ordered so that Sj(X) = Sj(X) 2 for all j = 1,2, ...,n. Let 
p : {1, 2, . . . , n} — > {1, 2, . . . , n} (usually not onto) satisfying the follow- 
ing properties: 

(a) Ej = S p(j) , 

11 



( b ) UJ=i s ; = Uj=i s p(j), and 

(c) E p (j) = Sp( fc ) if and only if p(j) = p(k). 

The map p is a choice of the indices of those, and only those, different Ey. 
Let c\ := max{p(j)}™ =1 and C2 := max[{p(j)}™ =1 \ {ci}]. We define the set 

fp:={AGC: s Cl (X) = s C2 (X)}. 

By virtue of Proposition Fp has empty interior. Furthermore, if all 
eigenvalues of -P(A) have geometric multiplicity equal to 1, then 

T P = {A G C : s n (A) = s n _!(A)}. 

Proposition 10 // all the eigenvalues of P(X) have geometric multiplicity 
equal to 1, then either Tp = or Tp lies on an algebraic curve (including 
the possibility of isolated points). 

Proof. Let 

T = {XeC: Sj (X) = s k (X), j^k} 

so that Tp C T . This set is the locus of all points (x, y) G M 2 such that 
the discriminant of the real polynomial in S defined by (jlljl is zero. The 
hypothesis ensures that J^C, and thus, either T = or T is an algebraic 
curve. The result follows just because Tp is a subset of T . □ 

In particular, Tp might include straight lines, single points, the empty 
set, or be a complicated set such as a Voronoi diagram (see Example 1 below). 

Borrowing a geological term, we call the set Tp the set of fault points of 
-P(A). In general, Tp will be made up of fault lines. The explicit determina- 
tion of the fault lines of P(X) requires computations with determinants and 
discriminants, and is therefore unrealistic. However, the following consider- 
ations demonstrate the role that the fault lines frequently play in the study 
of pseudospectra. 

Let F e (x,y) be as in (|12|). As mentioned above, apart from A = 0, if 
VF(i, y) does not exist, then x + iy G Tp. At these points, the curve 
dA £ (P) will typically fail to have a tangent line. There are other points 
where the tangent line will be undefined, those where \7F(x, y) = 0. In this 
case, there is a saddle point in the minimal singular value surface. These 
may or may not lie on Tp (see Section EJ). 
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Figure 1: The thin solid lines are Tp. The thick solid lines are dA e (P). 

Example 3 below illustrates a case in which Tp is a singleton. In Exam- 
ple 4, Tp is empty but there is, nevertheless, a point at which <9A £ (P) has 
no tangent. In Examples 1, 2, and 5, Tp is, indeed, made up of fault lines. 

Example 1 Let A be an n x n normal matrix with eigenvalues {Aj}™ =1 . 
Then the fault lines of P(X) = IX — A (i.e., for the standard eigenvalue 

the boundary 

□ 



problem) form the Voronoi diagram defined by {Aj}™ =1 (i.e. 
of their Dirichlet tessellation). 



Example 2 Naive experiments with diagonal matrix polynomials provide 
an insight on the possible structure of individual fault lines. For instance, 
let P(X) = diag{A 2 - 2A, (a - A)(A + 2)} and set w(x) = 1. In Figure HJ 
we depict the evolution of the set Tp and <9A £ (P) (e = 1/a/2, 1, a/3) for 
a = — 1, i, 1/2, 1. The fixed eigenvalues of -P(A) are plotted as "+" and the 
perturbed eigenvalue a is marked with a "o" . 

In general, an unbounded Tp appears to be more likely to occur. Nonethe- 
less this set can also contain a compact smooth curve. In both of the lower 
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Figure 2: dA E (P) has no tangent line at the origin for just one e. 

figures, Tp consists of an unbounded curve, which is asymptotic to a ver- 
tical line, and a closed compact curve on the right half plane around the 
perturbed eigenvalue a. The curve dAi(P) has a self intersection at A = 1 
for a = —1, i, 1/2. This can be shown from the fact that this part of the 
pseudospectrum depends only on the first diagonal entry of -P(A). The self 
intersection disappears as soon as a moves sufficiently close to 1. There are 
critical values of a, where 1 G Tp. Two of these critical values are a = 2/3 
and a = 4/3. □ 

Example 3 Tp can also be a singleton. In the left part of Figure |2l we 
depict dA e (P) for the linear matrix polynomial 



P(X) 



A + 3i/4 
A 





1 



-5/4 1 

A + 3/4 



the weight function w(x) = 1 and e 2 = 1/10, 5/16, 1/2. The very special 
structure of this matrix polynomial ensures that Tp = {0}. The boundary of 
the pseudospectrum does not have a tangent line at A = when e = a/5/16. 
Compare with Example 4 below. □ 

By construction, Tp is independent of w(x). Therefore, the singularities 
occurring on dA £ (P) in places where the gradient of (fT2J) fails to exist, are, 
with the possible exception of A = 0, independent of the chosen weights. In 
order to illustrate this remarkable fact, we consider two more examples. 

Example 4 The set Tp might be empty but the smoothness of dA e (P) 
might be broken at A = due to the weight function. Indeed, let n — 1, 
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y o- 




Figure 3: The thin solid lines are Tp. The thick solid lines are dA £ (P). 

P(A) = (A - l) 2 and w(x) = 2x + 1. Then si(x + iy) = (x - l) 2 + y 2 and 
T P = 0. 

When £ — 1, Fx(x,y) = x 2 -\- y 2 — 2x — 2^/x 2 + y 2 . Hence, F x (x,y) = 
if and only if 

(x - l) 2 + y 2 > 1 and y A + 2(x 2 - 2x - 2)y 2 + (x 4 - 4x 3 ) = 0. 

Thus, the curve cM.! (P) has a parameterisation of the form 



V±K X ) 



± J 2 + 2x - x 2 - 2V2x~TT ; -1/10 < x < 



in a neighbourhood of the origin. As d x y + (0) < and d x y-(0) > 0, G 
dAi(P) is a singularity of Lipschitz type. The boundaries of A e (P) for 
£ = 1/2, 1, 3/2, are drawn in the right part of Figure 121 □ 

Example 5 Let P(A) = diag{A 2 -l, A 2 -2A}. In FigureH we depict T P and 

<9A £ (P) for w(x) = 1 and e = a/3/5, 1, 2 (left), and for w(x) = x 2 + x + 1 
and e 2 = 1/20, 1/10, 1/5 (right). Here, Tp comprises a circle centred at 
(1/2,0) and the line x = 1/2. As in the previous examples, "+" marks the 
locations of the eigenvalues of P(X). □ 

All the above examples were designed in such a manner that both the fault 
points and the boundaries of pseudospectra can be constructed analytically 
either by hand or using algebraic computer packages. We produced FiguresHJ 
12] and El using commands provided in the standard distribution of Maple. 
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5 On the number of connected components 

Consider an n x n matrix polynomial -P(A) as in (Q), a real e > 0, and a 
weight function w(\) with w(0) = Wo > 0. Theorem 2.3 of [12] will be useful 
in the remainder of the paper. First we examine the case in which a(P) 
contains multiple eigenvalues more carefully and without the restriction of 
boundedness. A technical lemma will assist in the argument. 

Lemma 11 Suppose A and E are two nxn complex matrices such that the 
determinants detA and det(A + E) are nonzero. Then there is a continuous 
map t i-» E{t) G C nxn , t G [0, 1], such that E(0) = 0, E(l) = E, and 

det{A + E(i)) + and \\E(t)\\ < \\E\\ ; te [0,1]. 

Proof. Since detA ^ and det(A+E) ^ 0, no eigenvalue of the pencil A+tE 
can be equal to or 1 (and some may be infinite). So it may be assumed 
that det(v4 + tE) has s real zeros in the interval (0, 1), where < s < n. 

If s — 0, then the continuous map t \— > tE, t e [0, 1], has the properties 
required by the lemma. If s > 1, then let t± < ti < ■ ■ ■ < t s denote the zeros 
of det(A + tE) in (0, 1). For any tj (j = 1, 2, . . . , s), the matrix A + tjE 
is singular and for 8j > sufficiently small, we have 

det[A + (tj + e*%-)£] ^ and ||(^ + e 4 %-)£|l < ll^ll I G [0, 2tt]. 

In [0, 1], we replace each interval [tj — 8j,tj + 8j] with the circular arc 

C s = {t,-e%:0G[O,7r]}, 

and consider the continuous curve 

S= [0,t 1 -Sx]UC 1 U[t 1 +6i,t2-S 2 ]UC2U---U[t a - 1 +S a - U t s -S a ]UC a U[t a +6s,l] 

in the complex plane. For every continuous map t i— > z(t) G S, t G [0, 1], 
such that z(0) = and z(l) = 1, the map t h+ z(t)E G C nxn , t G [0,1], 
has the required properties. □ 

We are now ready to establish our main result on the number of con- 
nected components of pseudospectra. We should remark that, when A e (P) is 
bounded, the following theorem is a consequence of Theorem |21 Indeed, since 
w(x) is a real polynomial, u»(|A|) is a subharmonic function in C so, by The- 
orem|21 s(A) _1 iy(|A|) is subharmonic in C \ <j(P). If A £ (P) had a connected 
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component where there is no eigenvalue of P(A), then s n (A) w(|A|) _1 would 
have a local minimum in this component, which is impossible according to 
Theorem |21 

Theorem 12 If the matrix polynomial P(X) has exactly k (< nrri) distinct 
eigenvalues (not necessarily simple), then for any e > 0, the pseudospectrum 
A e (P) has at most k connected components. 

Proof. If A e (P) = C, then there is nothing to prove. So assume that 
A e (P) 7^ C, and consider a perturbation 

Q(X) = (P m + A m )A m + ■ • • + (P 1 + Ax)A + P + A 

in B(P,s,w) with det(P m + A m ) 7^ 0. By Lemma ITT1 there is a continuous 
map t h-> A m (t) G C nxn , t G [0, 1], such that A m (0) = 0, A m (l) = A m , and 

det(P m + A m (t)) 7^ and ||A m (t)|| < ||A m || ; tE [0,1]. 

Hence, every member of the family 

Q t (X) = (P m + A m (t))A m + ... + (P 1 + tA,)X + P + tA ; i G [0, 1] 

has exactly nm eigenvalues, counting multiplicities. Moreover, all Qt(X) 
(t G [0,1]) belong to B(P,e,w). Their eigenvalues lie in A e (P) and trace 
continuous curves from the eigenvalues of P(A) (= Qo(X)) to the eigenvalues 
of Q(X) (= Qi(X)). Thus, as in the proof of Theorem 2.3 of [12], the set 

A = {fiEC: det Q(ji) = 0, Q(X) G B(P, e, w), det(P m + A m ) ^ 0} 

has at most k connected components determined by the k distinct eigenvalues 
ofP(A). 

Now let Ao be an interior point of A e (P), and let R(X) = Y^T=o Rj ^ 
be a perturbation in B(P,e,w) with detP m = 0, such that Ao G a(R). 
Since A £ (P) 7^ C, R(X) has less than nm (finite) eigenvalues and, without 
loss of generality, we may assume that R(X) G dB(P, s n (A ) u>(|A |) -1 , w) C 
lnt[B(P, e, w)] (see Lemma |Hj). Then Ao is also an eigenvalue of all matrix 
polynomials 

R a (X) = (R rn + aI)X m + R m ^X m - l + --- + R 1 X + R -(aX™)I; a G C\{0}, 
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where det(P m + al) ^ and R a (X) G B(P,e,w) for sufficiently small \a\, 
i.e., A lies in A . By Corollary El A e (P) does not have more connected 
components than Ao (C A e (P)). Hence, A e (P) has at most k connected 
components. □ 

In this theorem, recall that since the leading coefficient of P(A) is non- 
singular, for e sufficiently small, A e (P) has exactly k bounded connected 
components. Thus, our upper bound for the number of connected compo- 
nents of A e (P) is always attainable when A e (P) is bounded. 

Proposition 13 If A e (P) is bounded, then any Q(\) G B(P,e,w) has an 
eigenvalue in each of these components. Furthermore, P(A) and Q(X) have 
the same number of eigenvalues (counting algebraic multiplicities) in each 
connected component of A £ (P). 

Proof. See Theorem 2.3 of HI. □ 



6 Multiple eigenvalues of perturbations 

In this section, we obtain necessary conditions for the existence of pertur- 
bations of P(A) with multiple eigenvalues. However, we first construct two 
perturbations of P(A) in B(P,e,w), which are of special interest. They are 
used in an argument generalising that of Alam and Bora (Theorem 4.1 of 
|2|) for the standard eigenvalue problem. 

Suppose that for a fi G A e (P) \ o~(P), the (nonzero) minimum singular 
value of the matrix P(/x) has multiplicity k > 1. Let also 

Sl(fi) > S 2 (/i) > •• > Sn-k(fj) > S n -k+l{M) = ••• = S„(/i) >0 

be the singular values of P(/x) with associated left singular vectors U\, u 2 , ■ ■ ■ , u r , 
and associated right singular vectors t>i,t>2, • • • ,v n . These singular vectors 
satisfy the relations P(p,)vj = Sj(fjL) Uj for j — 1, 2, . . . , n. 

Define the n x n unitary matrix Z = \u\ Ui • • • u n ] [v\ v 2 ■ ■ ■ v n ]* and 
the nxn matrix Z = [u n _ k+1 u n _ k+2 ■■■ u n ] [v n - k+1 v n „ k+2 ■■■ v n \* of rank 
k. Then Zvj = Uj for all j = 1,2, ... ,n, Zvj = Uj for all j = n — k + 
1, n — k + 2, . . . , n, and Zvj = for all j = 1, 2, . . . , n — k. Furthermore, 
the (nonsingular) matrix E = — s n (/j,) Z satisfies 

(P(h) + E)vj = s n (fi)uj — s n (fi)uj = 0; j = n — k+l,n — k+2, . . . ,n (14) 
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and 

(P(/i) + E)vj = sj(n) Uj - s n (n) Uj ^ ; j = 1, 2, . . . , n - k. (15) 

Similarly, the (rank k) matrix E = — s n (/x) Z satisfies 

(P(ji) + E)vj = ; j = n-k + l,n-k + 2,...,n (16) 

and 

{P{n)+E) Vj = Sj {n) Uj ^ 0; j = 1,2,..., n-k. (17) 

Note also that \\E\\ = \\E\\ = s„(u). 

M ii ii ii > L \r / 

Now define (for a given weight function w(x)) the matrices 



A ^ = (r-y^Hlvir'E; j = 0,1, 



m 



and 



7i XJ 



Aj = (-^) w.wdfiiy'E; j=0,l,...,m, 
I A* I 



where we set /x/ |/x| = when /x = 0. Then 

m / m \ 

^4V = ^^| Ai p'U(H)- 1 E = js 

i=o \i=o / 

and 

m / m \ 

^4V = l^2w s \i t \ j )w(\ ft \)- 1 E = E. 

3=0 \j=0 J 

Thus, for the (full rank) perturbation of P(X) 

Q(X) = (P m + A m )X m + ■ • • + (A + A X )A + P + A (18) 

and the (lower rank) perturbation of P(X) 

Q(X) = (P m + A m )X m + ... + (P 1 + Ai)A + P + A , (19) 

we have Q(p) = P{^) + E and Q(jjl) = P{fi) + E. From (UK), dHJ), ffifi and 
(|17p. it is clear that /i is an eigenvalue of the matrix polynomials Q(A) and 
Q(A) with geometric multiplicity exactly k and associated right eigenvectors 

v n-k+l, v n-k+2i ■ ■ ■ i v n- 
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Moreover, for every j = 0, 1, . . . , m, 



|Aj = wjwdfiir'WEW = » < e Wj 



Wj s n (jx) 
w(\n\) 



and 



|aj = Wjwdniy'WEW = w 7^ < eWi . 



w(\n\) 
Consequently, Q(A) and Q(A) lie in B(P,e,w) and the next result follows: 



Proposition 14 Let // G A e (P) \ cr{P) and let the nonzero singular value 
s n {n) (< £w i\f i \)) °f the matrix P(/J,) have multiplicity k > 1. Then the 
perturbation Q(X) in (|18|) and the perturbation Q(X) in (JT5)l lie in B(P,e,w) 
and have [i as an eigenvalue with geometric multiplicity equal to k. 

Clearly, every fault point of P(X) in C \ c(P) is a multiple eigenvalue of 
Q(X) and Q(X) with geometric multiplicity greater than 1. Furthermore, in 
the above discussion, note that for every j — n — k + 1, n — k + 2, . . . , n, 

u*P((S) = s n (fS)v* 

and 

u*{P{n)+E) = u*P{n)-s n {ii)u*Z 

= s n (fi) v* - s n (fj,)(Z*Uj)* 

= s n (fi) v* - s n (fi) v* = 0. 

Similarly, for every j — n — k+ 1, n — k+ 2, . . . , n, we have u* {P(f^) +E) = 0. 
Thus, u n -k+i, M„_fc+2, ■ ■ ■ ,u n are left eigenvectors of the perturbations Q(X) 
and Q(X) in (|TH|) and (fE?|) . corresponding to /i. 

The perturbations Q(X) and Q(X) defined by (fTH|) and (fE?|) depend on 
w(x) (which is considered fixed) and on the choice of [i. It is also worth 
noting that for /i = and a given weight function w(x) with a constant 
coefficient wo > 0, the construction of Q(X) and Q(X) is independent of the 
non-constant part of w(x) and requires only Wq. In the remainder of this 
paper, and without loss of generality, for the definition of Q(X) and Q(X), we 
use the constant weight function w c (x) = wq (> 0) instead of w(x) whenever 
// = 0. 

Using Lemmaim one can estimate the (spectral norm) distance from P(X) 
to the set of matrix polynomials that have a prescribed fi ^ o~(P) as an 
eigenvalue (cf. Lemma 3 of [To]). 
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Corollary 15 Suppose \x ^ cr(P), and let 5 = s n (//)iu(|//|) _1 . Then the 
perturbations Q(X) and Q(X) lie on dB(P, 5, w) and have \x as an eigenvalue. 
Moreover, for every e < 5, no perturbation of P(X) in B(P,e,w) has \x as 
an eigenvalue. 

Proposition 16 Let fj, G o~(P), and let u,v G C n be left and right eigen- 
vectors of P(X) corresponding to \i, respectively. If the derivative of P(X) 
satisfies u*P'(fi)v = 0, then ft is a multiple eigenvalue of P(A). 

Proof. If the geometric multiplicity of fi G o~(P) is greater than 1, then the 
proposition obviously holds. Hence, we assume that /i is an eigenvalue of P(A) 
with geometric multiplicity 1. For every vector y G C n , u*P(fi)y = 0, and 
thus, u _L Range[P(/x)]. Since u _L P'(n)v and the dimension of Range[P(//)] 
is n — 1, it follows that the vector P'(/i)v belongs to Range[P(/i)], i.e., there 
exists a y^ G C n such that 

P(fi)y, + P'(fi)v = 0. 

This shows that \x is a multiple eigenvalue of P(A) with the Jordan chain 
{v, y^} (see jH] for properties of Jordan chains of matrix polynomials). This 
implies that \i is a defective multiple eigenvalue of P(A). D 

Recall the function F e (x,y) = P e (x + iy) (x,y G R) defined in ()12|). 

Proposition 17 Suppose that for a point \x = x^ + iy^ o/ A e (P) \ o~(P), 
s n (fj,) is a simple singular value of P(/i) and u^v^ are associated left and 
right singular vectors, respectively, assuming that w(x) = w c (x) (= wo > 0) 
when /i = 0. Let 5 = s n (/i) w(|yu|) _1 (< e) and consider the perturbations 
Q(A), Q(X) G <9i3(P,o",u;) defined by flUD and (HHJ) . J/ £/je gradient of the 
function F$(x,y) = P^x + iw) ai /i zs zero, £nen ^ is a defective eigenvalue 
°f 0(A) and Q(A) TOi/j geometric multiplicity 1. 

Proof. Suppose /i / 0, and let VF^ij,,^) = 0, or equivalently (see 
Lemma EJ) , let 

Re ^^^^J =^^^ and Re [ U »^y- V ») =6 ^y-- 

Since 

-^- =P(a*) and _ =»P( A ,) I 
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we see that 



Moreover, 



*.«^0-H.(^ 



dWm " ?>'(H) and ^M = r>'OI), 



9a; |/i| ' ch/ |/i 

and consequently, 



Q^M t , _ ^ dw(H) ^ dw(|/z|) _ s [i 
dx M <9x dy \n\ 



m ;pVK = «;-r^^ = *^p-i*^p = *o«' , (N)- 



Consider the perturbation 

Q(A) = (P m + A m )A m + • • • + (A + Ai)A + Po + Ao 

in (|18|). Then Q(X) lies on the boundary of the (compact) set B(P,S,w) C 
£>(P, e, w) and its derivative satisfies 

(m 
, /I //i is / * A \ w'(|/i|) 7 1 

|/i| ** w(|/x|) \n\ 

m MM) w 

= S^-w'dfiD-S^w'dfil) = 0, 

where w M and v^ are left and right eigenvectors of Q(X) corresponding to [i, 
respectively (see Proposition (TJ] and the related discussion). The same is 
also true for the perturbation Q(X) in (JT§j) and its derivative. By Proposi- 
tions HU and CHI /i is a multiple eigenvalue of Q(X) and Q(X) with geometric 
multiplicity 1. 

For [i = 0, the proof is the same, keeping in mind that the constant 
weight function w c (x) = wo (> 0) is differentiable (with zero partial deriva- 
tives) at the origin. □ 
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7 Multiple points on dA e (P) and connected 
components of A £ (P) 



At first glance it may seem that multiple (crossing) points on dA e (P) will 
be exceptional. However, when we consider the evolution of dA £ (P) as e 
increases, it is clear that, as disjoint components of A e (P) expand, there will 
be critical values of e at which they meet and multiple points are created. 

Next, based on the results of the previous section, we show that multi- 
ple points of dA E (P) are multiple eigenvalues of perturbations of -P(A) on 
dB(P, e, w) and, also, these perturbations can be constructed explicitly. (Re- 
call that, when /i = 0, we use the constant weight function w c (x) = w > 
for the definition of the perturbations Q(X) and Q(X) in (fTHj) and (JUJ).) 

Theorem 18 Suppose that, as the parameter e > increases, two different 
connected components of A e (P) ^ C, Q\ and Q 2 , meet at \i G C Then the 
following hold: 

(i) If fj, t^ 0, then it is a multiple eigenvalue of the perturbations Q{X), Q(X) G 
dB(P,e,w) defined by flE) and (JTHJ). 

(ii) If /j, — and w(x) = w c (x) (= U7 > 0), then \i = is a multiple 
eigenvalue of the perturbations Q(X),Q(X) G dB(P,e,w c ). 

(iii) If fi = 0, w(x) ytz w c (x), A e (P) is bounded and the origin is the only 
intersection point of Q\ and Qi, then ft = is a multiple eigenvalue of 
a perturbation on dB(P, e, w) . 

Furthermore, in the first two cases, if s n ([i) is a simple singular value 
of P(n), then [i is a defective eigenvalue of Q{X) and Q(X) with geometric 
multiplicity 1. 

Proof. Suppose that s n (fi) (= ew(\fi\)) is a multiple singular value of 
the matrix P(fi). Then by Proposition H31 the perturbations Q(X),Q(X) G 
dB(P, e, w) have \i as a multiple eigenvalue of geometric multiplicity greater 
than 1. Hence, we may assume that s n (n) is a simple singular value of P(/i), 
and consider the three cases of the theorem. 

(i) Suppose /i 7^ 0, and recall fTHj) . By virtue of Lemma IH F e (x,y) is real 
analytic in a neighbourhood of /i. Furthermore, VP e (/i) = 0, otherwise the 
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implicit function theorem would ensure the existence of a smooth curve on 
a neighbourhood of \x parameterising dA e (P) and contradict the fact that 
dQ\ PI cX?2 is a finite set (Theorem |7J). Therefore, Proposition 1171 yields the 
desired conclusion. 

(ii) If /i = and w(x) = w c (x) (= w > 0), then the result follows by 
applying Proposition IT7I as in case (i). 

(iii) Suppose A e (P) is bounded, w(x) ^ w c (x), and /x = is the only 
intersection point of Q\ and Qi- By Proposition [T31 for any positive 5 < e, 
all the perturbations in B(P, 6, w) have a constant number of eigenvalues in 
A.g(P) C\Qj, say Kj, for j = 1, 2. Here and throughout this proof, eigenvalues 
are counted according to their algebraic multiplicities. 
Define the sets 

B = {Q(X) G B(P,e,w) : G a(Q)} C 9i3(P,e,w) 

and 

Bj = {Q(A) G i3 : Q(A) has less than Kj eigenvalues in Qj \ {0}} ; j = 1, 2. 

If Bj = (j = 1, 2), then PropositionUHland the continuity of the eigenvalues 
of matrix polynomials with respect to the entries of their coefficients imply 
that ^ Qj\ this is a contradiction. Hence, the sets B\ and B2 are both 
non-empty. 

Now consider the constant weight function w c (X) = w (> 0) and the 
associated e-pseudospectrum of -P(A), 

KvAP) = {/iGC: detg(/i) = 0, ||A || <e% Ai = --- = A m = 0}. 

Clearly, A £)Wc (P) C A e (P) and G 9A £)Wc (P). For any j = 1,2, consider a 
perturbation 

Qi(A) = (P m + A m )A m + --- + (P 1 + A 1 )A + P + A g S,-, 

and define the matrix polynomial 

<2i,c(A) = P m A m + • ■ ■ + Pi A + P + A = P(A) + A G B n 9z5(P, e, w c ) 

and the continuous trajectory 

Q,(£; A) = (P m + tA m )A m + ■ ■ ■ + (P ± + tA ± )X + P + A G B ; < t < 1 
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with Qj(0;X) = Qj fi {X) and Qj(l;X) = Qj(X). If // = is a multiple 
eigenvalue of Qj(t; A) for some t G [0, 1], then there is nothing to prove. 

Let fi = be a simple eigenvalue of Qj(t; A) G £> for all t G [0, 1]. Since 
A e (P) is bounded and the origin is the only intersection point of Q\ and Q 2 , 
by the continuity of the eigenvalues with respect to the coefficient matrices, 
it follows that all Qj(t;X) (0 < t < 1) have exactly Kj — 1 eigenvalues 
in Qj \ {0}, i.e., they lie in By Thus, Qj >c (X) G £>j. Moreover, again by 
Proposition EH and the continuity of eigenvalues, an eigenvalue of the matrix 
polynomials P(X) + (1 — t)A (0 < t < 1) traces a continuous path in Qj 
connecting the origin with an eigenvalue of -P(A). This means that the origin 
is an intersection point of A £<Wc (P) fl Q\ and A eu , c (P) fl Q 2 . Hence, // = is 
a multiple point of dA £>c (P), and as in (ii), it is a multiple eigenvalue of the 
perturbations Q{X), Q{X) G dB{P, e, w c ) C 9B(P, e,w). □ 

Now we can generalise a theorem of Mosier concerning scalar polynomials 
(Theorem 3 of [Tl]). 

Theorem 19 Suppose A e (P) zs bounded and Q is a connected component 
of A E (P). T/ien i/ie matrix polynomial P(X) has more than one eigenvalue 
in Q (counting multiplicities) if and only if there is a perturbation Q(X) G 
B(P,e,w) with a multiple eigenvalue in Q. 

Proof. For the converse part, it is clear that if a perturbation Q(X) G 
B(P,s,w) has a multiple eigenvalue in Q, then by Proposition IT3*1 P(A) has 
at least two eigenvalues in Q, counting multiplicities. 

For the sufficiency, if the matrix polynomial P(A) has a multiple eigen- 
value in Q, then there is nothing to prove. Thus, we assume that P(A) has 
two simple eigenvalues, Ai and A2, in Q. By the continuity of the eigenvalues 
with respect to the coefficient matrices, it follows that there is a positive 
S < e, such that A,y(P) has a (bounded) connected component Qg C Q that 
is composed of two compact sets, Qxj and Q 2 ,8, with disjoint interiors and 
intersecting boundaries. Moreover, without loss of generality, we can assume 
that Ai and A2 lie in the interior of Q± t s and Q 2 ,s, respectively. Then the 
curve enclosing Q$ either crosses itself or is tangent to itself at some point 
|ij6C\ °"(P)- The result follows from Theorem 1181 Note that if Q\^ and 
Q 2> s intersect at the origin and one other point, Theorem ITH1 does not apply 
at \x = 0, but it will at that other point. □ 
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8 Two numerical examples 



We present two numerical examples, which illustrate the results of the previ- 
ous section and suggest possible applications. The figures were drawn using 
the boundary-tracing algorithm described in |12j . 

Example 6 The spectrum of the 2x2 quadratic matrix polynomial 



P(A) 



(A -I) 2 




A 

(A-2) 2 





"10" 


■» o 


" -2 


1 




"10" 


= 




A 2 + 






A + 






1 







-4 




4 



is cr(P) = {1,2}. Both eigenvalues are plotted as "+" in Figure |H and have 
algebraic multiplicity equal to 2 and geometric multiplicity equal to 1. The 
boundaries dA e (P) for w(x) = x 2 + x + l (i.e., for perturbations measured in 
the absolute sense) and for e = 0.005, 0.0091, 0.02, 0.03, are also sketched 
in Figure EJ 

Assuming that the pseudospectrum A .oo9i(P) is connected with one self- 
intersection fj, = 1.4145 plotted as "o", this figure indicates that A e (P) 
consists of two connected components for e < 0.0091, and that it is connected 
for e > 0.0091. Moreover, the singular values of the matrix P(/x) are Si(/x) = 
1.4650 and s 2 (n) = 0.0402, i.e., s 2 (fi) is simple (/i is not a fault point of P(A)) 
and the function 

Po.oo9i(z,y) = F omgl {x + iy) = s 2 {x + iy) -0.0091 w(\x + iy\) ; x,y <ER 

has zero gradient at the point \x. Thus, by Proposition El and Theorem E3 
two perturbations of P(A) on the boundary of B(P, 0.0091, w) that have \x as 
a defective eigenvalue are Q(X) and Q(X) in (|18|) and (|19|) . and can be easily 
constructed. Left and right singular vectors of P{[J>) corresponding to si(/i) 
are 



Ml 



0.9726 
0.2325 



and V\ 



0.1141 
0.9935 



respectively, and left and right singular vectors of P(f^) corresponding to 
S2(/x) are 



u 2 = 

respectively. 

The unitary matrix 



-0.2325 
0.9726 



and t> 2 



Z = [ui u 2 ][vi V 2 ] 



-0.9935 
0.1141 



0.3419 0.9397 
-0.9397 0.3419 
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Figure 4: A single intersection point. 

satisfies 

Zv\ = Mi, u\Z = v\, Zvi = u 2 and u\Z = v%, 

and the rank one matrix 



u 2 v 2 



0.2310 -0.0265 
-0.9663 0.1110 



satisfies 

Zv\ = 0, u[Z = 0, Zvi = u<i and u\Z = v* 2 . 

We define the matrices 



A = Ai = A 2 = ( fi 2 + fi+ i)-\-s 2 ( f i)z) 

and the matrices 

A = A 1 = A 2 = ( fI 2 + fI + i)-\s 2 ( f i)z) 



-0.0031 -0.0086 
0.0086 -0.0031 



-0.0021 0.0002 
0.0088 -0.0010 



all with spectral norm 0.0091. Then the perturbations 

Q(X) = P(A) + (A 2 A 2 + A 1 A + A ; 
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0.9969 -0.0086 
0.0086 0.9969 



A 2 + 



-2.0031 
0.0086 



0.9914 
-4.0031 



A+ 



and 



<5(A) 



0.9979 0.0002 
0.0088 0.9990 



A 2 + 



P(A) + (A 2 A 2 + A 1 A + A ) 

A + 



-2.0021 1.0002 
0.0088 -4.0010 



0.9969 -0.0086 
0.0086 3.9969 



0.9979 0.0002 
0.0088 3.9990 



lie on dB(P, 0.0091, w) and have a multiple eigenvalue (approximately) equal 
to \i — 1.4145 with algebraic multiplicity 2 and geometric multiplicity 1, 
confirming our results. □ 

It is important to note that, by Theorems ITHl and \T§\ pseudospectra yield a 
visual approximation of the distance to multiple eigenvalues, i.e., the spectral 
norm distance from an n x n matrix polynomial -P(A) with a nonsingular 
leading coefficient and all its eigenvalues simple to n x n matrix polynomials 
with multiple eigenvalues. For a given weight function w(x), this distance is 
defined by 

r(P) := min{e > : 3 Q(\) G B(P,e,w) with multiple eigenvalues} 
= minje: > : 3 Q(X) G dB(P,e,w) with multiple eigenvalues}. 

Then Theorems E3 and E3 imply the following result (see [3 EI] for the 

standard eigenvalue problem). 

Corollary 20 Let P(A) be an n x n matrix polynomial as in (£Q) with a 
nonsingular leading coefficient and simple eigenvalues only. 

(a) If A e (P) is bounded, then 

t(P) = minje: > : A £ (P) has less than nm connected components} . 

(b) If A e (P) is unbounded and, as e increases from zero, its connected com- 

ponents meet at points different from the origin, then 

t(P) = min{e > : the number of connected components of A e (P) decreases}. 
Example 7 Consider the 3x3 self-adjoint matrix polynomial 



P(A) 



10 0" 




2 


A 2 + 


5 






3-1 
0-16 



A + 



-1 





3 








10 
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Figure 5: A damped vibrating system. 



(see [T2"l Example 5.2]), which corresponds to a damped vibrating system. 
The boundaries of A e (P) for w(x) = ||A 2 ||x 2 + ||Ai||x + \\A Q \\ = 5x 2 + 
6.3x + 10 (i.e., for perturbations measured in a relative sense) and for e = 
0.02, 0.05, 0.1, are drawn in Figure El The eigenvalues of -P(A), —0.08 ± 
il.45, -0.75±i0.86 and -0.51±il.25, are plotted as "+". 

We learn from this figure and the above discussion that there exist an 
£i = r(P) in (0.02, 0.05) (for which, the pseudospectrum starts having less 
than six connected components) and an e 2 in (0.05,0.1) (for which, the 
pseudospectrum becomes connected) such that the following hold: 



1. For every e < e\, all the perturbations Q(X) G B(P,s,w) 
simple eigenvalues. 



have only 



2. For every e G [£1,62), some perturbations Q(X) G B(P,e,w) have 
multiple non-real eigenvalues (in a neighbourhood between the eigen- 
values of -P(A) in the open upper half-plane and in a neighbourhood 
between the eigenvalues of -P(A) in the open lower half-plane), but no 
perturbation in B(P, e, w) has multiple real eigenvalues. 



3. For every e > e 2 , some perturbations Q(X) G B(P,e,w) 
real eigenvalues in the interval [—2.1, —0.2]. 



have multiple 

□ 
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